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ABSTRACT 

Wc present empirical relations between the global dust reddening and other 
physical galaxy properties including the Ha luminosity, Ha surface brightness, 
metallicity and axial ratio for star-forming disc galaxies. The study is based on 
a large sample of ~ 22 000 well-defined star-forming galaxies selected from the 
Sloan Digital Sky Survey (SDSS). The reddening parameterized by color excess 
E{B — V) is derived from the Balmer decrement. Besides the dependency of 
reddening on Ha luminosity / surface brightness and gas phase metallicity, it is 
also correlated with the galaxy inclination, in the sense that edge-on galaxies are 
more attenuated than face-on galaxies at a give intrinsic luminosity. In light of 
these correlations, we present the empirical formulae oi E{B — V) as a function of 
these galaxy properties, with a scatter of only 0.07 mag. The empirical relation 
can be reproduced if most dust attenuation to the H ii region is due to diffuse 
interstellar dust distributing in a disc thicker than that of Hii regions. The 
empirical formulae can be incorporated into semi-analytical models of galaxy 
formation and evolution to estimate the dust reddening and enable comparison 
with observations more practically. 

Key words: galaxies:ISM-galaxies:abundance~Hli regions-dust, extinction. 



1 INTRODUCTION 

Dust is a crucial component of galaxies in modifying the 
observed properties of galaxies by absorbing and scatter- 
ing starlight, and also re-emits the absorbed energy in 
the mid- and far-infrared bands. The extinctiorQ cross 
section of dust generally decreases with increasing wave- 
length, i.e., the extinction is more severe at shorter wave- 
length, particularly at ultra-violet (UV) band. As a re- 
sult, the spectral energy distributions (SED) of galaxies 
appear redder, which is so-called "reddening" . The effects 
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^ The nomenclature conventionally adopted is that, "extinc- 
tion" describes the absorption and scattering of light from a 
point source behind the dust screen; the "attenuation" de- 
scribes the reduction and reddening of the light in a galaxy 
where stars and dust are mixed in geometrical distribution. 
In this work, we do not discriminate between these two terms 
when referring to the dust attenuation within the galaxy, ex- 
cept for the Galactic "extinction curve" and the "attenuation 
curve" for the starburst galaxies. The extinction and redden- 
ing can be converted to each other with an assumed extinction 
curve. 



of dust should be properly accounted for when interpret- 
ing the observations of the galaxies, for example the lu- 
minosity or star formation rate (SFR) of a galaxy. For 
a galaxy with spectroscopic data covering suitable wave- 
lengths, it is usually possible to estimate the amount of 
dust extinction from the spectrum. For example, the av- 
erage dust reddening of H ll regions can be derived from 
Balmer decrement, and be used to correct for the dust at- 
tenuation to the emission-line luminosity. In other cases, 
one has to rely on empirical relations to get an estimate 
of dust extinction in a statistical way. 

Additionally, it is necessary to incorporate dust into 
the galaxy formation and evolution model. Currently 
there are two ways to account for dust effects: full radia- 
tive transfer calculations assuming dust properties and 
geometric distribution, or adopting simple recipes ob- 
tained empirically. Since the dust formation in the galac- 
tic environment is a complicated chemical process, it has 
not yet been implemented in the current galaxy evolution 
models. The dust properties and distribution are input 
in the models and should be tested with observations. 
Alternatively, many authors have adopted the empirical 
relations to compute the perpendicular optical depth of 
a galactic disc, and then assign a random inclination an- 
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gle for each galaxy to get the final dust corrections (e.g. 
Guiderdoni & Rocca-Volmerange 1987; Kauffmann et al. 
1999; Somerville & Primack 1999; De Lucia, Kauffmann 
& White 2004; De Lucia & Blaizot 2007; Kang et al. 2005; 
Kitzbichler & White 2007). Therefore, a well-defined em- 
pirical recipe of dust reddening will help the comparison 
of model predictions with observations. 

Previous studies have suggested that dust reddening 
in star-forming galaxies is correlated with the SFR, which 
can be well estimated with the intrinsic Ho? luminosity 
(e.g. Kennicutt 1998; Calzetti et al. 2010). Galaxies with 
large SFR (or high luminosity) show strong global ex- 
tinction in the emission lines (e.g., Wang & Heckman 
1996; Chariot & Fall 2000; Calzetti 2001; Stasinska & 
Sodre 2001; Afonso et al. 2003; Kewley et al. 2004; Zo- 
ran, Barkana, & Thompson 2006; Calzetti et al. 2007; 
Garn et al. 2010), and also large infrared to ultraviolet 
fiux ratios (e.g., Iglesias-Paramo et al. 2006 and reference 
therein) . 

In addition, it has been recognized that dust redden- 
ing might also be a function of gas-phase metallicity, i.e. 
the reddening and extinction increase with metallicity. 
For instance, the average reddening to H ll regions in an 
individual galaxy is correlated with the metallicity in the 
disc (e.g., Quillen & Yukita 2001; Boisser et al. 2004). 
For the case of whole galaxies, the relationship between 
reddening and metallicity also exists (e.g. Heckman et al. 
1998; Buat et al. 2002; Asari et al. 2007). It is also sup- 
ported by the fact that the low-metallicity galaxies such 
as blue compact galaxies are usually less reddened (e.g. 
Kong 2004). 

Recently, Garn & Best (2010) also found a signifi- 
cant correlation between dust extinction and metallicity, 
however, they claimed that the dependency of reddening 
on stellar mass is more fundamental. They built a sample 
of about 90 000 star-forming galaxies using Sloan Digital 
Sky Survey (SDSS) data, and compare the relationship 
between the dust extinction of Ha luminosity and SFR, 
metallicity as well as stellar mass, respectively. They con- 
cluded that the dust extinction Ahq can be best predicted 
from the stellar mass, with a scatter of ~0.3 mag. 

Besides, the dependency of dust attenuation on in- 
clination has been also investigated. While the correla- 
tion between dust attenuation of optical stellar contin- 
uum and disc inclination has been well established for 
disc galaxies (e.g. Driver et al. 2007; Shao et al. 2007; 
Unterborn & Ryden 2008), that of emission lines is still 
uncertain (e.g. Yip et al. 2010). 

These works focused on how the dust reddening de- 
pends on one property of galaxies, luminosity, metallicity 
or stellar mass, which are generally correlated with each 
other (e.g. Garnett & Shields 1987; Skillman et al. 1989; 
Zaritsky et al. 1994; Zahid et al. 2011; Tremonti er al. 
2004). Disentangling the dependencies of dust reddening 
on these galaxy properties can not only help us to derive a 
more accurate empirical formulae for dust reddening, but 
also provide more insight to the dust formation within 
galaxy environment. In this work, we study a large sam- 
ple of star-forming disc galaxies in the local universe with 
a median redshift of 0.07, selected from the spectroscopy 
database of SDSS. We obtain empirical formulae for dust 
reddening as a function of intrinsic Ha luminosity / sur- 



face brightness, gas-phase metallicity and the axial ratio 
of the disc. 

In Section[2]we describe our sample selection, and in 
Section[3]describe the spectral analysis and methods used 
to estimate the parameters such as the dust reddening, in- 
trinsic Ha luminosity, surface brightness, metallicity and 
axial ratio etc. We present the resultant expression for 
the empirical formulae of dust reddening, and compare 
our results with previous work in Section |4] In Section [5] 
we reproduce the observed trend with a toy model of 
parallel-slab disc, which gives some insights into the ge- 
ometry of dust distribution and dust-to-gas ratio. Con- 
clusions are given in Section |6] Throughout the paper, 
we will assume a cosmology of Ho = 72 km s~^Mpc~^, 
Q.M = 0.3 and Q.a = 0.7. 



2 DATA AND SAMPLE SELECTION 

We start from the spectroscopic sample of galaxies in the 
SDSS data release five (DR5, Adelman-McCarthy et al. 
2007). SDSS produces imaging and spectroscopic survey 
with a wide-field 2.5m telescope at Apache Point Ob- 
servatory, New Mexico (York et al. 2000). The survey 
provides imaging in five broad bands u, g, r, i, z, with 
magnitude limits of 22.2 in r band, and spectroscopic 
targets are selected using a variety of algorithms, includ- 
ing the "main" sample of galaxies with r-band magnitude 
brighter than 17.77 (Strauss et al. 2002) with fibers of 3'.' 
diameter. The spectra range from 3200 to 9200 A at a res- 
olution R = A/ A A — 1800 Q DR5 spectroscopic area cov- 
ers 5740 deg"'^, and there are ~675,000 spectra classified 
as galaxies. To select clean galaxies, we use the SDSS pho- 
tometric fiags to eliminate the targets that are only one 
part of a large galaxy or a part of merging galaxies. We 
list the sample-selection cuts and corresponding number 
of remaining objects within each subsample in Table [T] 
Then the galaxy spectra are corrected for the Galactic 
extinction using the dust extinction map (Schlegel et al. 
1998) with an extinction curve of Fitzpatrick (1999) with 
Rv = 3.1. 

We select emission-line galaxies where the Ha has 
been detected at high significance, i.e., S/N>20. To in- 
clude galaxies of high reddening, we use a looser criterion 
on H/3: S/N>10. The emission lines are measured with 
the procedures described in details in Section [3.11 In or- 
der to use the ratios of [O ill] /H/3 and [Nil] /Ha on line- 
ratio diagnostic diagrams (e.g. BPT diagrams; Baldwin, 
Phillips & Terlevich 1981) for spectral classification, we 
also impose the criteria [Olll] S/N>10 and [Nil] S/N>5. 
It should be noted that the requirement of lower limit 
in the S /N ratio for [N ll] may drop some low metallicity 
galaxies, and that for [O ill] may miss high metallicity 
galaxies. We will justify the S/N criteria in Section [4.21 

Galaxies with significant broad Ha component are 
rejected. The broad-line Active Galactic Nuclei (AGNs) 
are defined as objects for which adding an additional 
broad component of Ha to the emission-line model can 
significantly improve the fit to the Ha+[Nll] blend (refer 
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to Dong et al. 2005, 2007; also Zhou et al. 2006). In prac- 
tice, the galaxies with broad Ha component detected at 
the ^ 5a significance level are regarded as candidates of 
broad-line AGNs, and removed from the sample. We also 
remove the narrow-line active galaxies by using the BPT 
diagram (Kauffmann et al. 2003b, hereafter Ka03; Kew- 
ley et al. 2006) based on [Nii]/Hq. The galaxies below 
the Ka03 pure star-formation line on [Nn]/Ha diagram 
are referred as our sample of star-forming galaxies. Most 
of these galaxies lie below the extreme-starburst line on 
the [S ll]/Ha and [O l]/Ha diagrams, with a fraction of 99 
percent and 94 percent, respectively. Using a more strict 
criterion on selecting Hll galaxies given by Stasiriska et 
al. (2006) will result in less metal-rich galaxies (see S I3.3p . 

As shown by Kewley et al. (2005), if the nuclear spec- 
trum contains less than twenty percent of the total galaxy 
light, we will likely over-estimate the global metallicity 
and reddening, and under-estimate the global SFR by a 
significant fraction. Therefore, we remove the galaxies for 
which the fiber magnitude at g-band is greater than the 
total magnitude by 1.7 mag. 

In the following analysis, we will pick out the late- 
type (presumably disc-dominated) galaxies from the star- 
forming galaxies based on the likelihoods provided by 
SDSS pipeline. The photometry pipeline provides the 
likelihoods (dev_L, expJL, and starJL) associates with the 
de Vaucouleurs, exponential, and PSF fits, respectively. 
The fractional likelihoods for the exponential fit is calcu- 
lated as 

exp_L 



/(exp_L) 



exp_L -I- dev_L -I- star_L ' 



(1) 



and similarly for /(exp_L) and /(star_L). It is suggested 
that the fractional likelihood greater than 0.5 for any of 
the three model fits is generally good as a threshold for 
object classification (Stoughton et al. 2002). For a galaxy, 
/(starJL) is generally zero. We define a galaxy as a disc 
galaxy if the logarithmic likelihood for an exponential 
fit (r-band) is larger than that of a de Vaucouleurs fit 
by 0.2 dex (corresponding to /(exp_L)> 0.6), while as 
an elliptical galaxy if the logarithmic likelihood for a de 
Vaucouleurs fit is larger than that of an exponential fit 
by 0.2 dex (/(dev_L)> 0.6). This parameter is also corre- 
lated with the compactness index such as R50 /Rgo, which 
is commonly used in quantitative classification (e.g. Shi- 
masaku et al. 2001; Strateva et al. 2001). With this cri- 
terion, 23919 star-forming galaxies are classified as disc 
galaxies, and 7650 as elliptical galaxies, and 595 as un- 
classified type. 



3 METHODS AND SAMPLE PROPERTIES 

3.1 SPECTRA ANALYSIS: 

STARLIGHT-CONTINUUM 
SUBTRACTION AND EMISSION-LINE 
FITTING 

In order to measure the emission lines, we take two steps 
to analyze the spectra: continuum fitting and emission- 
line fitting. First, we subtract the stellar continuum fol- 
lowing the recipe described by Lu et al. (2006). In brief. 
Ensemble Learning for Independent Component Analysis 



(EL-ICA) has been applied to the simple stellar popula- 
tion library (Bruzual & Chariot 2003, BC03) to derive a 
set of templates, which then are shifted and broadened 
to match the stellar velocity dispersion of the galaxy, and 
reddened assuming a starburst-like extinction law to fit 
the observed galaxy spectra. During the continuum fit- 
ting the bad pixels fiagged out by SDSS pipeline as well 
as the emission-line regions are masked. From the fit, 
we obtain simultaneously the modeled stellar-light com- 
ponent, stellar velocity dispersion and an effective red- 
deninflfl to the stellar light (refer to Lu et al. 2006 for 
details). Second, the emission lines are modeled with var- 
ious Gaussians on the continuum-subtracted spectra, us- 
ing the MPFIT package (Markwardt 20090 implemented 
in Interactive Data Language (IDL) . The formal la errors 
in fiux obtained from the fitting, are propagated from the 
error of the spectra, and then adopted as the emission- 
lines flux uncertainties. The emission lines we measure 
have been corrected for absorption lines, by subtract- 
ing the stellar component models. The typical absorption 
correction is 27 percent of the flux of H/? emission-line. 
We examine the model-fitting of higher-order Balmer ab- 
sorption lines, like H(5, which is less contaminated by 
emission lines, and estimate the uncertainty in absorp- 
tion measurement to be less than 24 percent, including 
the statistical uncertainty. Thus the uncertainty in H/3 
absorption measurement is generally less than 6 percent 
of H/? emission flux. 

The corresponding emission-line regions are masked 
in the continuum fits. To determine the proper mask- 
ranges for the emission lines, usually several itera- 
tions of the above procedures are required. Emission 
lines, Hq, H/3, H7, US, He, [On]A3727 A, [Neiii] 
AA3869,3969 AA, [Nev] A3427 A, [Oiii]A4363 A, [Hen] 
A4686 A, [O III] AA4959, 5007 AA, [Ni] A5199 A, [Hei] 
A5876 A, [Nii]AA6548,6584 AA, [S ii]AA6716, 6731 AA, 
[Oi]AA6300,6344 AA, [Arm] A7136 A are included in 
the first iteration, but insignificant ones (signal-to-noise 
ratio in emission-line fiux, S/N<3) are dropped in the 
later fitting. For robustness, we assume identical profiles 
for [Nil] doublet lines and Ha. The [Sll] doublet lines 
are assumed to have a same profile, so are [O ill] doublet 
lines. The ratios of [N 11] doublets and [O ill] doublets are 
fixed to their theoretical values, 2.96 and 3, respectively. 
[0 11] doublets are each modeled with a single gaussian of 
the same width. To reduce the uncertainty in the mea- 
surements of weak emission lines, we fix their profiles to 
those of strong lines of similar ionization states. As the 
final procedure, upper limits are given to the undetected 
lines assuming that the lines have the identical profile as 
the detected strong lines. If no emission-line has been de- 
tected significantly (> 5a), no emission-line fiux will be 
given for that spectrum. 



^ The effective reddening derived in this process is fairly well 
correlated with that of emission lines estimated using the 
Balmer decrements for Hll galaxies (see also e.g. Calzetti et 
al. 1994; Stasiriska et al. 2004). 

* MPFIT package includes routines to perform non-linear least 
squares curve fitting, kindly provided by Craig B. Markwardt, 
available at jhttp:/ /purl.com/net/mpfit| 
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Table 1. Summary of sample-selection cuts. 



Selection cut Num remained Percent removed by cut 



SDSS DR5 spectroscopic sample 


582 


512 


/ 


Photometric clean sample 


495 


165 


1 A nn 

14. yy 




191 


289 


61.37 


S/N(H/3) > 10 


128 


046 


33.06 


S/N([Oiii]) > 10 


67 


236 


47.49 


S/N([Nii]) > 5 


67 


013 


0.33 


Remove broad-line AGN candidates 


57 


068 


14.84 


Remove multiple observations 


56 


241 


1.45 


Select star-forming galaxies using BPT 


46 


865 


16.67 


Require >20% total light in the fiber 


32 


164 


31.37 


Select disk galaxies 


23 


919 


25.63 


Require -2.5 < N2 < -0.4 


22 


616 


5.45 



3.2 REDDENING AND Lh^ CORRECTION 

We estimate the reddening of emission lines using the 
Balmer decrement Ha/H/3 ratio. An intrinsic value of 
2.86 as expected for case-B recombination with electron 
density — lOOcm"'' at T = 10* K is assumed (Os- 
terbrock & Ferland 2006). This value is generally consis- 
tent with the lower limit of the measured Hq/H/3 ratio 
(Figure [TJi,) in our sample. There is only a small fraction 
(about 1.5 percent) of objects with Ha/H/3 below 2.86, 
likely due to measurement uncertainty. For these objects, 
the reddening is adopted as zero. Due to our stringent cri- 
teria for H/3 and [O in] detections, a significant fraction 
of objects with high Balmer decrement values have been 
dropped; thus most of our sample have Ha/H/3 < 7. The 
attenuation A^a is estimated from the Balmer decrement 
Ha/H/3 assuming an extinction curve, and then used to 
correct Ha luminosity. 

We also apply aperture correction on Ha luminos- 
ity, based on the difference between the model magni- 
tude and the fiber magnitude at p-band. This correction 
method assumes the distribution of Ha emission is the 
same as that of the stellar hght (continuum emission). 
Such a correction is only an approximation because the 
line emission and the stellar continuum emission may not 
distribute in the same way. Some other authors use em- 
pirical approach to make aperture corrections taking into 
account the color differences within / outside the fiber 
(Brinchmann et al. 2004), or constrain the global SFR 
from fitting stochastic models to the photometric SED 
(Salim et al. 2007). The former method is based on the 
main assumption that the distribution of specific SFR 
for a given set of colours inside the fibre is similar to 
that outside. We test with this method, but find that at 
a given set of colours the likelihood distribution of spe- 
cific SFR inside the fiber varies with Balmer decrement. 
The typical specific SFR is higher for galaxies with larger 
Balmer decrement, and lower for galaxies with smaller 
Balmer decrement. Thus this method of aperture cor- 
rection may introduce dependence of Ha luminosity on 
Balmer decrement. Since our goal is to investigate the 
correlation between dust reddening and luminosity, we 
decide to settle for the simple scaling method. We also 
examine if the global SFR obtained with the method of 
Salim et al. (2007) is used, and compare with the SFR 



estimated from far-infrared luminosity, as we will check 
for our corrected SFR in the following part of this sec- 
tion. The test suggests our simple method is no worse 
than theirs. 

Because different parts of a galaxy suffer from dif- 
ferent extinction, the extinction derived from the Balmer 
decrement is only a certain average. In the following, we 
will check if the aperture and attenuation correction in- 
troduce any fake correlation between the Ha luminosity 
and reddening. If the corrected /(Ha) does not represent 
an accurate intrinsic Ha fiux, then the star-formation 
rate estimated from Ha luminosity will be inaccurate. To 
examine this issue, we use far-infrared luminosity as the 
reference tracer for SFR and compare it with the SFR 
estimated from Ha luminosity. On one hand, although 
integrated IR emission should provide a robust measure- 
ment of SFR in dusty circumstances (Kennicutt 1998 and 
references therein. Dale & Helou 2002), there are also 
calibrations based on luminosities of specific bands at in- 
frared (e.g. Wu et al. 2005; Alonso-Herrero et al. 2006; 
Zhu et al. 2008; Calzetti et al. 2007, 2010; Rieke et al. 
2009). At high luminosity (Lro^m > 1.4 x 10''^ ergs s"^), 
L-jOf_im correlates linearly with SFR, thus could be used 
as a tracer of SFR (Calzetti et al. 2010, their equation 
22): 

SFR{70) (Mgyr"') = 5.88 x W'^'^Lro^.m- (2) 

On the other hand, we adopt the calibration of Calzetti 
et al. (2010, their equation 5) to convert Ha luminosity 
to SFR as: 

SFR{Ra) (Mgyr"') = 5.45 x 10'*^LHa(ergs s~'), (3) 

in which Lna should be corrected for intrinsic extinction. 
This calibration is based on solar metallicity and Kroupa 
(2001) Initial Mass Function (IMF). The Kroupa IMF 
has two power laws, one with a slope of —1.3 for stellar 
masses ranging from 0.1 to 0.5 Mq and the other with 
a slope of —2.3 for stellar masses ranging from 0.5 to 
100 Mq. This calibration is based on a t > 1 Gyr age 
constant star-formation stellar population. 

We cross-match our parent sample of star-forming 
galaxies with the 70/im band photometry catalogs 
from Spitzer Wide-area InfraRed Extragalactic Survey 
(SWIRE; Lonsdale et al. 2003) Data Release 3. To esti- 
mate the Ha SFR for the galaxies with confidence, we 



© RAS, MNRAS 000,[T]-?? 



Dust reddening in SF galaxies 5 




Figure 1. Distributions of properties for the selected ~22 000 disc-dominated H II galaxies in the SDSS DR5 (see text): (a) Balmer 
Decrement Ha/H/3; (b) dust-attenuation corrected Ho luminosity log i^SI' ("^) 12 + log(0/H); (d) axial ratio h/a; (e) Ha surface 
brightness corrected for dust attenuation log , see Equation (6] (f) stellar mass. 




Figure 2. The ratio of SFR{Ha)/SFR(70) versus Balmer decrement Ho/H/3 in logarithmic scale. SFR(H«) (MQyr^l) = 5.45 X 
10^'*'^LHa(ergs s~^), where Lhq is attenuation-corrected with a starburst attenuation curve (Calzetti 2001, left panel), or a 
Galactic extinction curve (middle panel: O'Donnell 1994 curve, right panel: Fitzpatrick 1999 curve). SFR(70) (A/Qyr^^) = 
5.88 X 10~**L70(jm. Solid lines represent the best-fitting for the sample, and dashed lines represent SFR(Ha)/SFR(70) = 1. 



require the aperture to include at least twenty percent of 
the total light (Kewley et al. 2005). We also require the 
detection of Ha emission to be more significant than 10a. 
In order to get a matched sample of reasonable size, we 
apply a looser S/N criteria on other emission lines (i.e. 
S/N>10 for Ha, and S/N>5 for [Nil], H^, [Olll]). With 
a matching radius of 5", we get 156 star- forming galaxies 
with measurements of 70/xm fluxes. Removing five galax- 
ies with contaminating sources nearby in IR emission, two 
galaxies with unreliable or problematic 70/im fluxes, and 
one galaxy with problematic aperture correction (neg- 
ative value), there are 148 galaxies left in the SWIRE 
star- forming galaxy sample. 

We correct Lua for dust attenuation with three ex- 
tinction curves respectively: 1) the attenuation curve for 
the continuum of starburst galaxies (Calzetti et al. 2000); 
and Galactic extinction curve of 2) O'Donnell (1994); or 
3) Fitzpatrick (1999). Converting Lhq to SFR with Equa- 
tion |3l we investigate the ratio of SFR(Ha)/SFR(70) as 



a function of the Balmer decrement in the log-space. If 
Lua is properly corrected for dust attenuation, the ra- 
tio of LHc</L70Mm, thus the ratio of SFR(Hq)/SFR(70) 
should be independent of the dust reddening and the 
Balmer decrement. In the SWIRE star- forming galaxy 
sample, there are 147 galaxies with Lropm greater than 
1.4 X lO"*^ ergs s"^ ~ 3.7 x 10** Lq, for which Equa- 
tion [2] can be used to estimate SFR(70). Figure [2] shows 
that if the starburst attenuation-law is adopted (left 
panel), SFR(Hq)/SFR(70) still correlates with Ha/H/S 
(the spearman rank coefficient Vs = 0.33, with the prob- 
ability of null hypothesis Pnuii < 10~^), implying that 
the Lhq might have been over-corrected than demanded, 
while the Galactic extinction-curves give better correc- 
tion on average (middle and right panel). With Fitz- 
patrick (1999) curve, the ratio SFR{Ua) / S F R{70) is un- 
correlated to the Balmer decrement (r^ = —0.01, Pnuii = 
0.90). This means the H« SFR can be well determined 
from the attenuation-corrected Lua by assuming Fitz- 
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Patrick's curve. Therefore, we will adopt this curve for 
attenuation correction to Lhu in the following analysis. 

For our final sample of star- forming disc galaxies, the 
corrected Lhq is shown in distribution histogram in Fig- 
ure [TJb). The dust-extinction corrected Ha luminosity is 
in the range from 4x lO^'^ ergs s~^to 2x 10*^ ergs s~^, with 
a median value of 3x10*^ ergs s^^. The typical error in 
Z/Ha is 0.04 dex, as the quadrature sum of the measure- 
ment uncertainty in observed flux of Ha emission, and 
the Balmer decrement used for attenuation correction. 
The uncertainty induced by either the average reddening 
we simply assumed or the aperture correction for Lhq 
has not been included. Assuming Fitzpatrick's extinction 
curve, the color excess is estimated as 

E{B -V) ^ 1.99 log "^/gf^ for Ha/H/3 ^ 2.86 

= for Ha/H/3 < 2.86. (4) 

The median formal uncertainty in Ha/H/3 for our se- 
lected sample is typically 3.4 percent, which gives an un- 
certainty of 0.03 mag in E[B — V). The uncertainty of 
E{B — V) is about 0.04 mag typically for the small sam- 
ple of SWIRE star- forming galaxies. We also calculate 
E{B-V) from the ratio of Ha/H ^ assuming Fitzpatrick's 
extinction curve, and the derived E{B — V) is quite con- 
sistent with the value obtained using Equation 2] This 
proves that our corrections for Balmer absorptions are 
quite robust. 

We note that in the relation between 
SFR(Ha)/SFR(70) and the Balmer decrement, ir- 
respective of which extinction curve adopted, there is 
a moderate scatter ~0.13 dex. This scatter includes 
the measurement errors of observed Lua, the Balmer 
decrement and model/fiber magnitudes at g-band; the 
uncertainty in aperture-correction and attenuation- 
correction; the measurement error in observed I/7o,jm; 
and the calibration error in the SFR(Ha)/SFR(70) 
ratio. The overall measurement error for SFR(Ha) is 
typically 0.04 dex, and for SFR(70) typically 0.01 dex. 
Then the remaining scatter of 0.12 dex accounts for 
the sum in quadrature of uncertainties in aperture and 
attenuation correction, and the the calibration scatter 
in SFR(Ha)/SFR(70). Therefore, any one of these 
uncertainties, for example the the calibration scatter in 
SFR(Ha)/SFR(70) should be less than 0.12 dex, which 
is smaller than the calibration uncertainty in SFR(70) 
of Equation |2l --0.2 dex (Calzetti et al. 2010, refer to 
their Section 5). Note in passing, we obtained a larger 
scatter (~0.17dex) in SFR(Ha)/SFR(70) if MPA-JHlf| 
SFR(Ha) is used instead. 

3.3 METALLICITY 

Estimating global metallicity has been widely studied by 
using strong lines. However, there is still no consensus on 
which line ratio should be used. Ratios of strong lines, 
such as R23 = log([Oii]A3727 A + [O iii]AA4959, 5007 
AA) / H/3, [Nii]A6584 A / [Oii]A3727 A, N2 = 
log([Nll]A6584 A / Ha), and 03N2 = log([0 lll]A5007 A 

^ MPA-JHU catalogue and the total SFRs are provided on 
the website |http: / / www.mpa-garching.mpg.de/SDSS /DR7/ 1 



/ H/3) / ([N ll]A6584 A / Ha) are used to estimate the oxy- 
gen abundance (e.g., McGaugh 1991; Storchi-Bergmann 
et al. 1994; Denicolo et al. 2002; Kewley & Dopita 2002; 
Pilyugin 2003; Pilyugin & Thuan 2005; Pettini & Pagel 
2004; Tremonti et al. 2004; Liang et al. 2006; Shi et al. 
2006; Nagao et al. 2006). The commonly used metallic- 
ity indicator R23 is not used here, because it includes 
[0 11] emission line which is much prone to dust extinc- 
tion. Note that one purpose of our work is to investigate 
the relation between dust reddening and the metallic- 
ity, we should avoid any spurious correlation potentially 
caused by the reddening correction. The estimator A'^2 
and 03N2 are both not sensitive to the reddening cor- 
rection. But the measurement error of 03N2 (typically 
0.03 dex) is larger than that of A'^2 (typically 0.01 dex). 
Therefore A'^2 is preferred as our metallicity diagnostic. 

The relation of line ratio versus abundance can be 
generally calibrated using two different approaches: the 
empirical method relying on the electron temperature 
Te as a surrogate for metallicity (cooling increases with 
metallicity), or a comparison with photoionization mod- 
els. Note that different methods lead to systematic dif- 
ferences in the calibrated relations up to a factor of three 
or more in some extreme cases, apart from the limited 
range for the validation of the relations (refer to Kewley 
et al. 2008 for a detailed discussion). 

We estimate the metal abundance with N2 index of 
Pettini & Pagel (2004, hereafter PP04), further revised 
by Nagao et al. (2006) and Liang et al. (2006). PP04 
derived their formula from Hll regions with well deter- 
mined oxygen abundance based on Te-method, while Na- 
gao et al. and Liang et al. calibrated gas metallicity in 
gala:xies obtained using other methods such as Te-method 
or Bayesian-technique by comparing with the theoretical 
models of Tremonti et al. (2004). The three calibrations 
yield very similar estimates of [O/H] at low abundances 
(e.g. 12 -I- log (O/H) < 8.5), but they deviates consider- 
ably at high abundances (e.g. 12+log (O/H) > 8.5): PP04 
gives a lower abundance value than those of Liang et al. 
(2006) and Nagao et al. (2006). Since PP04 calibration 
does not extend to high abundance (most of their objects 
have 7.5 < 12 4- log (O/H) < 8.5), we adopt the calibra- 
tion of Nagao et al. (2006), expressed as 

^■2 = 96.641 - 39.941 xy + 5.2227 x j/^ - 0.22040 x , (5) 

where y = 12 + log (O/H). Equation[5]is valid within the 
metallicity range 7.15 ^ 12 + log(0/H) ^ 9.16, corre- 
sponding to -2.5 ^ N2 ^ -O.4I3 We estimate metal- 
licity for our sample with 7V2 for consistency, thus we 
reject those with N2 > —0.4 (only about five percent 
in our sample). Note that this selection has little effect 
on the results we obtain in this paper. The final sample 
consists of 22616 galaxies, with the oxygen abundance 
in the range of 7.64 < 12 4- log (O/H) < 9.16 and with 
a median value < 12 + log (O/H) >= 8.84 (Figure [l};) . 



^ Nagao et al. (2006) did not specify the valid range of N2 
for their calibration, we adopt the range of their observational 
sample used for calibration. We are prudent to not use the 
extrapolation when > —0.4, which saturates and yields 
quite high and unreasonable mctallicitics. 
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Figure 3. The reddening of emission lines parameterized by E{B — V) s.s a function of (a) attenuation-corrected Ha luminosity 
Lua (in units of ergs s~^), (b) oxygen abundance for our final sample of 22616 star- forming disc- dominated galaxies. The red line 
is the median value, while the blue lines enclose 68 percent (Icr) objects at a given log or 12 -|- log (O/H). (c) Distribution of 

E{B-V). 



The oxygen abundances could be converted to metallic- 
ity in units of solar metallicities (^0=0.02), adopting a 
value of 12 + log(0/H)Q = 8.66 (Asplund et al. 2004), 
and the range is 0.1 < Z/Zq < 3.2, with a median of 
< Z/Zq >= 1.5. 

The error in the metallicity includes measurement 
uncertainty and calibration uncertainty. For our sample, 
the metallicity is estimated using N2, and the typical 
measurement uncertainty in A'^2 is 2.5 percent, which 
corresponds only up to 0.02 dex in 12 + log (O/H). The 
main uncertainty in metallicity estimation comes from 
the scatter of the calibration of N2 indicator. Nagao et al. 
did not explicitly provide the scatter of their calibration 
(Equation [5| . However, PP04 proposed a linear calibra- 
tion of N2 and indicated the Icr uncertainty in Z is 0.18 
dex. The A'^2 calibration in Equation[5]differs significantly 
{AZ ^ 0.2 dex) from PP04's calibration only at metallic- 
ities 12-hlog (O/H) < 7.5 and 12-hlog (O/H) > 8.5. More- 
over, Kewley & Ellison (2008) compared various metal- 
licity calibrations, and concluded that metallicities es- 
timated from strong-line methods should be consistent 
within 0.15 dex. Therefore, we estimate the upper limit 
of the uncertainty in metallicity to be AZ ~ 0.15 dex. 

3.4 SURFACE BRIGHTNESS Bhc, DISC 
INCLINATION AND OTHER 
PROPERTIES 

The surface brightness of a galaxy is the flux received 
from a unit solid angle as it appears on the sky, then 
we define the intrinsic surface brightness Bua as the Ha 
luminosity per area within the half-light radius of the 
galaxy, i.e. 

Bhc = ^ (ergs s^^kpc"^), (6) 

where Rso is the Petrosian half-light physical radius (in 
kpc) in r-band, which can be obtained by the product 
of Petrosian half-light radius (in arcsec) provided by the 
SDSS pipeline, and the angular distance. The apparent 



surface brightness is dimmer than Bua with redshift by 
a factor of (1 + z)"^. The distribution of Bua is shown 
on (Figure [T^), and the formal uncertainty for Bua is 
about 0.05 dex. We will use this definition of Bua to com- 
pare with a toy model in next section. We use the axial 
ratio in the exponential fit of the galaxy {ABexp given 
by SDSS pipeline) as a surrogate for the disc inclination 
(Figure [T)i). While the axial ratio b/a close to 1 means 
that the disc galaxy is face-on, b/a close to means the 
disc is edge-on. The typical error of b/a is about 0.02. 
Note that Petrosian radius is defined as the radius at 
which the ratio (equals to some specified value, 0.2 for 
SDSS) of the averaged surface-brightness in a local an- 
nulus to the mean surface-brightness within it (Blanton 
et al. 2001; Yasuda et al. 2001). Thus the Petrosian ra- 
dius is not very sensitive to the disc inclination. However, 
the Ha surface brightness Bho estimated in Equation [6] 
is only a crude average approximation, because Hll re- 
gions likely occupy only a portion of the galaxy, and the 
distribution may be very inhomogeneous. 

The stellar mass M, (Figure [if) is calculated from 
template fits to the SDSS five-band photometry with the 
sdss_kcorrect package (Blanton & Roweis 2007) in IDL. 
The magnitudes are corrected to the AB system (-0.036, 
0.012, 0.010, 0.028 and 0.040 for the ugriz filters, respec- 
tively) and corrected for Galactic extinction. In short, 
Blanton & Roweis (2007) built their five templates using 
a technique of nonnegative matrix factorization based on 
a set of 450 instantaneous burst stellar population models 
(Bruzual & Chariot 2003) and 35 emission templates of 
MAPPINGS-HI models (Kewley et al. 2001). Specifically, 
the stellar population models include all six metallicities 
{Z from 0.0001 to 0.05) and 25 ages (from 1 Myr to 13.75 
Gyr) based on Chabrier (2003) stellar initial mass func- 
tion and the Padova 1994 isochrones. For each stellar pop- 
ulation model, three different dust models are assumed: 
no dust extinction, rv = 3 with Milky Way-like extinc- 
tion, or TV = 3 with Small Maggellanic Cloud-like ex- 
tinction. They assumed a homogeneous dust distribution 
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and shell geometry for the latter two dust models. The 
authors compared their measurements of stellar masses 
for SDSS galaxies with those obtained by Kauffmann et 
al. (2003a), and found that the two set of masses are 
consistent with each other, with a scatter of only 0.1 dex 
(refer to their Figure 17 in Blanton & Roweis 2007). 

In addition to these parameters, we also look into 
other properties based on the spectra. The 4000 A dis- 
continuity is due to the opacity from ionized metals, and 
its amplitude indicates the stellar population ages. In hot 
stars, the metal elements are multiply ionized and the 
opacity decreases, thus the 4000 A break strength is small 
for young stellar populations and large for old, metal- 
rich populations (Bruzual 1983; Balogh et al. 1999). The 
break discontinuity is generally defined as the ratio of 
the continuum red- and blue-wards 4000 A. We adopt 
the narrow definition D„(4000) introduced by Balogh et 
al. (1999) using bands 3850-3950 A and 4000-4100 A. 
D„(4000) is measured from the reddening-corrected mod- 
eled spectrum, which is obtained in the continuum fits 
described in Section TS. 11 instead of directly from the ob- 
served spectrum. An extensive test shows that the former 
has two advantages over the latter. One is that for the ob- 
served spectrum of a low S/N ratio, the model-fitting pro- 
vides a filter to the noises on the observed spectrum. The 
other is, on the modeled spectrum, the bias in D„ (4000) 
introduced by reddening can be corrected, about an offset 
of 0.03 larger for the typical reddening {E{B-V) = 0.23) 
of our sample. YL5 absorption line equivalent width (EW), 
which reaches a peak for the stellar population of age 
~ 0.1-1 Gyr (the Lick index 115^, Worthey & Otta- 
viani 1997; Kauffmann et al. 2003a), is measured from 
the modeled continuum as well. The electron density is 
estimated from the ratio of [S ll] AA6716, 6731 AA (Oster- 
brock & Ferland 2006). 



4 RELATION BETWEEN INTRINSIC 
REDDENING AND OTHER GALAXY 
PROPERTIES 

In this section, we will investigate the correlations be- 
tween the intrinsic dust-reddening and other observable 
properties or those properties that could be deduced from 
observations. The properties include the gas metallic- 
ity, the attenuation-corrected Hq luminosity or surface 
brightness, stellar mass, and the inclination of the disc. 

For the ~22 000 star-forming disc-dominated galax- 
ies in our final sample, the distributions of emission- 
line reddening, which is parameterized by E{B — V), axe 
shown as a function of Lhq (Figure[3^) or 12 -I- log (O/H) 
(FigureOs). On Figure|3}3, the metallicity estimated from 
N2 has a cut-off at 12 -|- log (O/H) = 9.16, which corre- 
sponds to N2 = -0.4. The histogram of E{B - V) dis- 
tribution is shown on Figure [3Ic). While the E{B — V) 
range is wide (0.0~1.0 mag), 95 percent of our sample is 
in the range of 0.05 to 0.6 mag. 

Figure [3] shows that the reddening increases both 
with increasing Lua and 12 + log (O/H). The intrinsic 
dust-reddening is clearly correlated with the Ha lumi- 
nosity, with a Spearman rank coefficient = 0.59 (the 
probability of null hypothesis PnuU < 10"^) for our sam- 
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Figure 4. Left panel: log Lhq versus E{B — V). is in 

units of ergs s~^. Black dots represent for the galaxies with 
axial ratio b/a < 0.9, and blue dots for the face-on galaxies 
{b/a > 0.9). Red lines are the best-fitting power-law for face- 
on galaxies in each metallicity bin. Right panel: The E(B — V) 
ratio versus the axial ratio b/a for each metallicity bin. The 
E(B — V) ratio is the ratio of E{B — V) in other galaxies to that 
deduced from the linear function of Lhci for face-on galaxies. 
The solid line presents the best-fitting. The three metallicity 
bins arc 0.5Zq, Zq and 2Zq, as marked on the left panel. 



pie of star- forming disc galaxies (Figure [SK). This re- 
lation has been known for a long time (e.g. Wang & 
Heckman 1996; Calzetti et al. 2001; Hopkins et al. 2001; 
Afonso et al. 2003). The reddening correlates with Hq 
surface brightness as well (not shown in the figure), with 
rs = 0.51 (Pnuii < 10^''). On the other hand, the redden- 
ing is also correlated with the gas-phase metallicity indi- 
cated by oxygen abundance 12 -I- log (O/H) {rs = 0.67, 
fnuu < 10~^, Figure Et), which has also been known 
previously (e.g. Heckman et al. 1998; Boissier et al. 2004; 
Asari et al. 2007). At low metallicity, the reddening in- 
creases slowly with [O /H] , but at high metallicity, it rises 
more sharply. The correlation between E{B — V) and 
metallcity is the strongest among all of the correlations, 
therefore we divide our sample into bins of metallicity 
to investigate the relations of reddening with other pa- 
rameters (Ha luminosity /surface brightness, axial ratio, 
stellar mass, etc.) at a given metallicity bin in the rest 
part of this section. 

We divide the sample into continuous metallicity 
bins. The bin is so chosen that a bin-size is no less than 
0.05 in 12 -I- log (O/H), with at least 400 objects in each 
bin. In order to isolate the inclination-dependence, we 
check the relation between E{B — V) and Lua for nearly 
face-on galaxies, i.e. (b/a) > 0.90, and then compare it 
with the rest of the galaxies. Clearly, face-on galaxies 
have systematically small E{B — V) values than the oth- 
ers at a given intrinsic Lua (Figure U left panel). We fit 
the relation of reddening with luminosity for these face- 
on galaxies using a power-law function. For a galaxy, the 
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E{B — V) ratio is defined as the ratio of E{B — V) of 
this galaxy, to the expected E{B — V) based on its Ha 
luminosity from the best-fitting power-law function for 
the face-on galaxies. The E{B — V) ratio is found to 
correlate with the axial ratio, which can be well fitted 
by {b/ay. Figure |4] shows logLna versus E{B — V) as 
well as and the best-fitting curve for the face-on objects 
on the left panel, and the dependency of the E{B — V) 
ratio on the axial ratio b/a on the right panel. These 
results are only illustrated for three metallicity bins: 
12 + log(0/H) ~ 8.36, 8.64 and 8.94 (i.e. ~0.5Zq, Zq 
and 2Zq, adopting [12 + log (0/H)]q = 8.66 ± 0.05, As- 
plund et al. 2004). The correlation between the E{B — V) 
ratio and the axial ratio becomes stronger as the metall- 
city increases, for example rs = —0.49 (Pnuii < 10~^) for 
the bottom metallicity bin 12 -I- log (O/H) ~ 8.94. We fit 
the data in each bin with a joint function in the form of 



E{B~V) =pLt,^^,,{b/ay 



(7) 



where Lua,42 = Lua/iW^^ ergs s~^). For each metallicity 
bin, the fitting parameters are obtained by minimizing 
in the fitting implemented by MPFIT, and the results are 
listed in Table [2] Only uncertainties in E{B — V) have 
been considered in the fitting, thus the errors of the fitting 
parameters are underestimated. 

Figure [5] shows p, q and 7 as a function of the metal- 
licity. We can see that p increases with the metallic- 
ity, and can be fitted with an exponential function of 
12 + log (O/H), i.e. a power-law function of metallicity 
Z, in the form of 



= 0.072 exp( 



0.7i±o.oi 



) 



= 0.18^ 



Z 
Z^ 



(8) 



where a: = 12 -f log (O/H) — 8. The power-law index q 
varies slightly between 0.08 and 0.17, and the average 
value is about 0.15. We therefore adopt the parameter q 
as a constant: 



(g) =0.15. 



(9) 



I7I first increases from about 0.1 with metallicity and 
then reaches to a constant (I7I ~ 0.4) at about 12 + 
log (O/H) ~ 8.5. We use a step function to fit 7 as a 
function of 12 + log (O/H), which is expressed as 



7 = -0.1 for 0.1 <x < 0.5 
= -0.4 for 0.5 s; a; < 1.2 



(10) 



Substitute p, q and 7 into Equation [T] we will get the 
overall function of luminosity, metallicity and inclination, 
and we call this {Lhc Z, b/a) model hereafter. To ver- 
ify the reliability of our empirical formulae, we define 
the E{B - V) residual as AE{B - V) = E{B - V)obs - 
E{B — V)m.odei, and show its distribution in Figure |S] 
(the first column, upper panel). The E{B — V) distribu- 
tion can be well fitted with a Gaussian profile, with width 
o — 0.068 mag. If we adopt a more stringent S/N criteria 
(S/N> 100) on Ha emission, the result will yield a sample 
of ~5000 star-forming disc galaxies. We repeat the pro- 
cedure above and find the scatter of AE{B — V) remains 
similar, a = 0.066 mag, indicating that the scatter is in- 
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Figure 7. Left panel: log -Bhq versus E{B — V). -Bua is in 
units of ergs s~^ kpc Right panel: The E(B — V) ratio 
versus the axial ratio b/a for each metallicity bin. The E{B — 
V) ratio is the ratio of E{B — V) in other galaxies to that 
deduced from the linear function of -Bjjci for face-on galaxies. 
The three metallicity bins are 0.5Zq, Zq and 2Zq, as marked 
on the left panel. The symbols and the lines represent the same 
as in Figure |4] 



dependent of Ha S/N, in the S/N range we consider (we 
will discuss in more details and show that our S/N cri- 
teria is reasonable in the following Section [4. 2|l . Besides, 
we find that AE{B — V) is uncorrelated with the param- 
eters including Lua , metallicity or axis ratio (Figure [SJ 
the right three columns), which suggests that the (Lhq, 
Z, b/a) model considered above reproduces the intrinsic 
reddening well in the disc galaxies. 

Note that at low metallicities the scatter in the rela- 
tion between E{B — V) ratio and b/a is larger than that 
at high metallicities (Figure [5] the right column). The 
power-law index 7 is close to zero at low metallcity, i.e. 
the relationship between E{B — V) ratio and b/a is al- 
most flat; while 7 ^ —0.4 when the metallicity increases 
to solar or super-solar value. This suggests that dust red- 
dening is not so sensitive to the axis ratio at low metal- 
licity as at high metallcity. It may be attributed to two 
possible reasons. One is that the relation at low metallic- 
ity might remain as strong as at high metallcity, but the 
limited dynamic- range of E{B — V) and the measurement 
uncertainties prevent us from recovering the relationship. 
The E{B — V) dynamic-range is only about 0.2 mag at 
low metallicity, and increases to 0.8 mag at high metallic- 
ity. The other possible reason is that the relation indeed 
becomes weaker at low metallicity. At present we can not 
distinguish which one is true, especially when the num- 
ber of face-on galaxies at low metallicity is limited. If we 
adopt 7 to be —0.4 for all metallicities, the distribution 
of AE{B - V) slightly changes. 

We do similar analysis on how E{B — V) varies as a 
function of Ha surface brightness and inclination in the 
metallicity bins, because the Ha (SFR) surface density is 
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(Lua.Z.b/a) model {B„^,Z,b/a) model 




8.2 8.6 9.0 8.2 8.6 9.0 

12 + log [0/H] 12 + log [0/H] 



Figure 5. Fitting parameters as a function of metallicity. (a) p, (b) q and (c) 7 are from Equation [3 and (d) p, (e) q and (f) 7 
are from Equation 1111 Those symbols without error bars indicates the error is smaller than the symbol size. In each panel, the 
best-fitting function is shown in blue solid line. 



Table 2. Coefficients of empirical fits to reddening in a function as in Equation [T] the 
empirical relation of reddening parameterized by color excess E(B — y), as a function of 
-Lhq a-iid the axial ratio 6/a in different metallicity bins. 
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believed to physically link with the gas surface density, 
which is closely coupled with dust. Besides, the gas sur- 
face density is provided by the current models of galaxy 
formation. Therefore, the E{B — V) as a function of Ha 
surface brightness will enable direct application of the 
empirical relation on the intrinsic reddening in the galaxy 
evolution models. Figure [7| shows Bua versus E(B — V) 
and the best-fitting curve of a power-law function for the 
face-on objects on the left panel, while the distribution 



of E{B — V) ratio against the axial ratio for the three 
selected metallicity bins on the right panel. In each met- 
allcity bin, we fit the data with a function in the form 
of 

EiB^V) =pBS„,,o(6/ar, (11) 

where BHa,40 = -BHa/(10"*'' ergs s^^kpc"^). All the fit- 
ting parameters are listed in Table[31 We also fit p with an 
exponential function of 12-|-log (O/H), q with a constant. 
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Figure 6. The first column shows the distribution of E{B ~ V) residuals in our sample for the two empirical models: the (I/hqi 
Z, b/a) model and the (Bhc, Z, b/a) model. /\E{B - V) is defined as /\E{B -V) = E{B - V)obs - E{B - V)model- N is the 
fraction of the total sample. The best-fitting Gaussian is over-plotted with the centre shown in dotted line for each model. The 
second to forth column shows AE(B — V) as a function of Lhq or ^Hq (in units of log ergs s~^ or ergs s' -ikpc-2), 12-hlog(0/H) 
and the axial ratio b/a, respectively. In each panel, the sample is binned in corresponding parameters, and the median values with 
Icr dispersions in each bin are shown. The gray dashed line represents AE{B — V) = 0. 



and 7 as a step function (Figure [5|. 
p = 0.092 '^^P( o,6i±o.oi ) 



= 0.27=^"" 

(q) = 0.23, 
and 

7 = -0.2 for 0.1 < a- < 0.5 

= -0.4 for 0.5 ^ a; < 1.2 



(12) 
(13) 

(14) 

With the model expressed in Equation llll hereafter 
(Bhci, Z, b/a) model, we calculate AE{B — V) and 
show the distribution in Figure |6l The distribution can 
be also fitted well with a Gaussian profile, with width 
a = 0.068 mag. AE{B — V) is uncorrelated with the 
parameters including Bua, the metallicity or the axial 
ratio, which suggests that the (Bua, Z, b/a) model can 
reproduce the intrinsic reddening well in the disc galaxies 
too. 

In order to investigate the origin of the scatter in the 
empirical relations, we examine the correlation between 
AE{B — V) and other galaxy parameters, such as stellar 
mass, the 4000 A break index, EW of Balmer absorption 
lines, the electron density in Hll region or concentration 
index. The AE{B - V) for either (Lhc, Z, b/a) model 
or (Bna, Z, b/a) model is not correlated with any of the 
parameters above, implying that the dust reddening in a 
disc galaxy is primarily determined by only three physical 
properties: luminosity (SFR) / surface brightness (SFR 
surface density), metallicity and disc inclination. 

We also attempt another multi-parameter model 
(M«, Z, b/a) with similar procedure described above, and 
find that the scatter of AE{B — V) is 0.072 mag, which is 
slightly larger than the two models above (Equation[7|and 
lll[) . However, AE{B — V) for this model is still weakly 
correlated with Lua (ts = 0.15, Pnuii < 10~^). This sug- 
gests that (Lho or SFR, Z) combination works better 



than (M, , Z) combination. Physically, at a given metal- 
licity, the amount of dust is proportional to the amount 
of cold gas, which is closely correlated with SFR. 



4.1 Comparison with previous work 

We compare our models with previous empirical rela- 
tions. Garn & Best (2010) investigate the dependency 
of dust extinction at Ho? on SFR, metallicity and stellar 
mass respectively, in star-forming galaxies based on SDSS 
database. They conclude that the stellar mass is the most 
fundamental parameter in determining the dust extinc- 
tion in the local Universe, and their relation between stel- 
lar mass and dust extinction seems to hold up to redshift 
of 1.5 (Sobral et al. 2011). They utilize different models 
to fit the extinction Aua, with each model taking into 
consideration one of the three parameters: SFR, M* and 
metallicity. They find that the M, model gives the best 
fit, with the minimum of scatter (0.28 mag) in the extinc- 
tion residual A Ana = A^^, - A^"'^"'-, less than the scat- 
ter for SFR model and metallicity model, both of which 
are 0.33 mag. When calculating the Aua, they adopt 
the Calzetti et al. (2000) dust attenuation-curve. We re- 
peat their analysis for our selected sample of disc galaxies 
and compare the scatter of AAua among different mod- 
els. In each model a low-order (third order) polynomial 
is adopted to fit the extinction, and the distribution of 
AAua is well fitted by a Gaussian, with widths a =0.27, 
0.28 and 0.28 mag for the M,, Lhc (SFR) and metallic- 
ity model respectively. Higher order polynomial fits will 
not decrease the scatter. We check the scatter of AAhc 
for our (I/Ha, Z, b/a) or {Bua, Z, b/a) model assum- 
ing Calzetti's attenuation law, and the Gaussian width is 
only a = 0.22 mag (corresponding to a = 0.065 mag for 
the AE{B - V) distribution). We also check for the (M«, 
Z, b/a) model, and find that the scatter of AtIhq is 0.24 
mag, which is better than the single-parameter models, 
but worse than the other two multi-parameter models. 
Our empirical relations predict dust reddening better 
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Table 3. CoefBcients of empirical fits to reddening in a function as in Equation 1111 the 
empirical relation of E(B — V) as a, function of Bhq and the axial ratio b/a in different 
metallicity bins. 
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8.46 


539 


0.223±0'OO6 
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-0.149±0-022 
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< 
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h log (O/H) 
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8.51 
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0.205±0"05 


0.229±O-008 
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5.34 


8.51 


< 


12 ~ 


h log (O/H) 
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8.56 
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0.242±0'006 
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8.61 
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< 
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< 
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6.58 


8.86 


< 
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h log (O/H) 
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8.91 


2456 


0,355±o.oo3 


0.189±0003 


-0.406±0-006 


6.77 


8.91 


< 


12 - 


h log (O/H) 


< 


8.96 


2551 


0.413±ooo3 


0194±o.002 


-0.339±0-005 
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12 - 


h log (O/H) 
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9.01 


2191 


0.479±o.oo3 
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_0.346±o.oo6 
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h log (O/H) 
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9.11 


1052 


0,567±o.oo4 


0.286±0003 
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12 - 


h log (O/H) 


> 


9.11 
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526±o.oo5 


0.244±0-004 


-0.399±0-009 


11.22 



than the single-parameter models, probably because we 
include disc inclination in the analysis. In order to investi- 
gate what dominates the relationships of the variables, we 
apply principle component analysis on the physical vari- 
ables. The derived principle components (PCs) are linear 
combinations of the variables, and represent for the direc- 
tions of maximum variance in the data. For Hq luminos- 
ity and surface brightness, E{B — V), M, and metalhc- 
ity, the quantities are normalized to zero mean and stan- 
dard deviation of unity, while for the axial ratio original 
value is used. We test with different subsets of variables 
adopted in each multi-parameter model as well as all of 
the variables, and find similar results. Take the case of 
(Lhc«, Z, b/a) model for example, the first two PCs con- 
tribute to about 80 percent of the total variance. In PCI, 
the variables Lua, Z and E{B — V) weight almost equiv- 
alently. This is consistent with what Garn & Best (2010) 
found. But in PC2 the axial ratio weights the most (-0.98) 
and E{B — V) weights the second (0.30). The result is con- 
sistent with what we previously found. Although the axial 
ratio only weakly correlated with E{B — V) {va = —0.13, 
^null < 10"^), after removing the E{B — V) dependence 
on Ha luminosity and metallicity, the negative correla- 
tion between b/a and E(B — V) emerges. 



4.2 Error analysis 

We have shown that the global reddening of the emission 
lines from a disc galaxy can be reasonably determined 
by its Ha luminosity or surface brightness, gas metallic- 
ity and disc inclination. The la scatter of AE{B — V) is 
about 0.07 mag, which is contributed by several sources, 
including the measurement uncertainty in the Balmer 
decrement, the error in the metallicity, the uncertainties 
in other variables and the intrinsic scatter in the relation. 
We will discuss the uncertainties in the following. 

Let us take a review of the uncertainties in the phys- 



ical quantities. The formal error for Lhq measurement is 
0.04 dex, and for B-na is 0.05 dex. As discussed in Sec- 
tion [321 the uncertainty induced by aperture-correction 
and attenuation-correction is less than 0.12 dex. The typ- 
ical error of bja is 0.02. The uncertainty in metallicity 
calibration dominates the uncertainty of metallicity, esti- 
mated to be ^ 0.15 dex. We use Monte Carlo simulation 
to estimate the error of E{B — V) for Equation [7] and 
111! taking the uncertainties for the best-fit parameters 
and the measurement errors for the physical variables. 
The errors of E(B — V) predicted by Equation [71 and [TTI 
are both about 0.056 mag. If taking the uncertainty in 
aperture-correction and attenuation-correction (0.12 dex 
at most) into account in Lhq and Bhq error, the simu- 
lated E(B — V) errors slightly changes to 0.058 mag. Re- 
member that the measurement error of E{B — V) is 0.03 
mag, the scatter of AE{B — V) (0.07 mag) can be almost 
explained by the uncertainties of the physical quantities. 
Considering all these uncertainties, E(B — V) predicted 
from the empirical relations are quite consistent with the 
observed values. 

From Figure [51 we can see that the AE{B — V) dis- 
tribution can be well described by a Gaussian profile, 
indicative of a random distribution. However, the intrin- 
sic scatter of the empirical relation is still unclear, be- 
cause of the large uncertainty in the metallicity calibra- 
tion. Therefore, to estimate the intrinsic scatter of the 
relation, it will be of help to decrease the uncertainty 
in the metallicity calibration and obtain the metallicity 
more precisely. One may find there is slight excess over 
the positive tail of the gaussian distribution. We examine 
the objects (about one percent of the whole sample) with 
AE{B -V) > 0.2, i.e. the observed E{B - V) is under- 
estimated with our empirical models by 0.2 mag, and find 
that these objects are prone to large Balmer decrement, 
low H/3 S/N and low H/3 EW. 

As previously mentioned, we require S/N> 20 for 
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H«, and S/N> 10 for H/? and [Olll] when selecting our 
sample in Section [J] The S/N criteria for Ha will ex- 
clude some galaxies with low SFR and low Ha EW. Thus 
the number of galaxies with low Ha luminosity and low 
SFR will be reduced. In fact, when the criteria S/N> 10 
is applied for H/3, most of the remaining sample {'^99.9 
percent) will have Ha S/N> 20. In order to investigate 
whether the H/3 and [Olll] S/N criteria bias our results, 
we attempt a looser S/N(> 5) criteria for both lines, 
and the sample is expanded by 50 percent to have 34815 
galaxies. Then we fit the data with both the (Lua, Z, 
b/a) model and {Bua, Z, b/a) model. The results show 
that the fit parameters are quite similar as those we have 
obtained in Section 2] but the scatter in AE{B — V) 
increases to about 0.08 mag. Note that 10 percent un- 
certainty in Ha/H/3 will result in an uncertainty of 0.086 
mag in E{B — V), which explains the increased scatter. 
Thus including the low S/N objects merely introduces 
more scatter in the relation. 



5 A TOY MODEL ON DUST EXTINCTION 
TO Hii REGIONS IN A DISC GALAXY 

Two different sources contribute to the extinction for H ll 
regions: dust associated with the individual Hll region 
and diffuse foreground dust in the galaxy. Let us first 
consider the dust in the precursor molecular clouds (MC). 
The intrinsic extinction in individual H ll region depends 
on the relative distribution of young massive stars, gas 
and dust in the star-formation regions, which are also 
closely related to the age and size of the precursor MC, 
as well as the stochastic formation of high mass stars for 
a given initial mass function. An H ll region is embedded 
in the cold gas/dust envelope before the cold gas/dust is 
blown away. The time scale for such a process is expected 
to depend on the size and density of the precursor cloud: 
longer for more massive and dense H ll region. Panuzzo et 
al. (2003) discussed the attenuation to the Ha emission 
line for different evaporating time scale, in comparison 
with the lifetime of massive stars that produce the ion- 
izing continuum. They simply adopt a central ionizing 
star in a spherical cloud, which corresponds to the maxi- 
mum attenuation case. If the time scale is short enough, 
most Hii regions are produced outside of the MC, and 
the main attenuation is due to foreground dust. In the 
case of long evaporation time, an H ll region is buried in 
the MC during most of their life time and the extinc- 
tion associated with the Hll region will be important. In 
a galaxy, all these cases are likely present. Furthermore, 
the young ionizing stars are not necessarily forming in the 
cloud centre. These complications make any prediction of 
the reddening in an individual cloud rather difficult. 

However, as far as a large number of H ll regions are 
concerned, the statistical behavior may be quite well de- 
fined. The Ha luminosity function of H ll regions in spiral 
galaxies can be described by double power-law functions, 
with a fiatter power-law below the break luminosity of 
38.6 (in units of ergs s~^, logarithmic), and a steeper one 
above it (e.g. Kennicutt et al. 1989; Bradley et al. 2006). 
In our sample, most galaxies have Ha luminosities above 
40 (in units of ergs , logarithmic) with a median of 



41.5, and the maximum is 43.3. Therefore, a galaxy con- 
sists of numerous (order of hundreds to millions of) H ll 
regions. The Hll regions follow certain distributions of 
ages and initial masses, determined by the star-formation 
processes, because the star-formation time in a galaxy is 
likely much longer than the lifetime of an individual H ll 
region, i.e. the life time of OB stars. If these distributions 
are the same for each galaxy and there are a large number 
of Hii regions, one may expect that the average extinc- 
tion within H ll regions should be similar in each galaxy. 
In this respect, it is not expected that the dust reddening 
depends on the Ha luminosity or the inclination of the 
disc. In contrast, for example, if the MC mass function 
and stellar IMF depend on the global properties of the 
galaxies, such as the gas surface density, then one may 
expect that the average extinction is correlated with the 
global properties. However, this is poorly known so far. 

Next, we consider the extinction due to diffuse inter- 
stellar dust in a galaxy. It is expected that this extinction 
depends on the gas surface density, inclination, metallic- 
ity of the interstellar medium and size of the galaxy. To 
simplify the model, we will adopt smooth distributions 
for the dust, gas and star- formation regions with the fol- 
lowing assumptions: 

(1) dust-to-gas density ratio is proportional to the gas 
metallicity (Issa et al. 1990), i.e., gas and dust are cou- 
pled: 

Pd oc pg X Z, (15) 

where pd and pg are dust density and gas density, respec- 
tively; 

(2) both gas and dust distributions in the vertical (z) 
and radial (R) direction can be described by an exponen- 
tial law (e.g. de Jong 1996; Thomas et al. 2004; Misiriotis 
et al. 2006), i.e., 

pg = Ps{0) exp{-z/hg) exp(-7?/i?g), (16) 

where Pg(0) is the gas density at z = and R — 0, hg and 
-Rg are the characteristic lengths in the vertical and radial 
directions for gas distribution. But the dust distribution 
may have different length scales due to radial metallicity 
gradient: 

Pd = Pd(0)exp(-z//id)exp(-i?/i?d), (17) 

where Pd(0) is the dust density ai z — and R = 0, hd 
and _Rd are the characteristic lengths in the vertical and 
radial directions for dust distribution; 

(3) the Hii distribution is also exponential in both ra- 
dial and vertical direction, but with different character- 
istic lengths (e.g. de Jong 1996): 

PHa = PHa(O) exp(-2//lHa) exp(-i?/_RHQ). (18) 

By integrating pHa over z direction, we yield a surface 
density of Ha, Shq- Specifically, the Ha surface density 
at i? = is 

Shq,o = / pua{0)exp(-z/hHa)dz + / pHaiO)exp{z/hiia)dz 

Jo J-oo 

= 2pHc{0)hnc (19) 

The gas disc is rather thin in general. As a result, 
dust attenuation can be modeled as plane-parallel except 
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at extremely high inclinations. Young stars are formed 
in dense molecular cores, which might be distributed in 
a plane thinner than the gas disc. For a plane-parallel 
model, the optical depth is expressed as 



exponential, we will get i?g = 1.4i?HQ. With above rela- 
tions and the conversion between Lna and SFR (Equa- 
tion [3|, we can write: 



Pdkyz/ cos i, 



(20) 



where fc^ is the opacity coefficient, and i is the inclination 

o 

angle, and i = represents the line of sight perpendicu- 

o 

lar to the galaxy disc, while i = 90 represents the line of 
sight parallel to the galaxy disc. The optical depth to Ha 
emission from the middle plane of the disc to the upper 
disc surface at 7? = is defined as 



THct.O ~ 



Pd{0)K exp{~h/hd)dh = pd{0)k„hd- (21) 



and the optical depth from that point to the lower disc 
surface can be expressed as 



T2{z) = 2THa,0 - Tl{z) 



(23) 



Light from the Hll regions at the position {z, R) 
needs to travel through the disc, from both upper and 
lower sides. With these descriptions, we can calculate the 
Ha luminosity. Ha surface brightness and average attenu- 
ation at Ha and H/3 by integrating the radiation transfer 
function over both vertical and radial directions within 
the disc. Then the Ha luminosity is 

Wp[-^ exp(-i?/i?d)] + exp[-^ exp(-i?/J?dJp) 

Substitue Equation 1181 1221 1231 y = exp(— z/h,Ha) and 
a = R/Rhc into Equation 1241 then we get 

Lua = TvR^^'EB_a,o dy adaexp{-a){ 



exp[ 
exp[- 



COS I 



+ 



exp{-a/rj)y 
exp(-a/r,)(2-yi/«)]}. 



(25) 



where ^ = -RHa,max/i?HQ is set by the threshold surface 
density of star formation iiHa.max = -Rg ln(Eg,o/Sth); 
and the threshold gas surface density is adopted as Eth — 
5Mq pc~^(Martin & Kennicutt 2001); Eg,o is the surface 
density of gas at _R = 0; 77 = Rd/Rua and k = hd/hua- 
In order to reduce the number of free parameters (EHafi, 
THctfi, cosi, rj and k), we use the empirical relation 
between the SFR surface density Espr and gas surface 
density Eg (known as Schmidt-Kennicutt law; Kennicutt 
et al. 1998): 

EsFR (MQyr"^kpc"^) 

= (2.5 ±0.7) xlO^"* 



1.4±0.15 

^ I ■ (26) 



, A-fQpc 

Reminding that the gas and H 11 distributions are both 



(27) 



It is symmetric between the middle plane to the upper 
and lower disc surface, thus the overall optical depth 
through the whole disk vertically at i? = should be 
2thc«,o. Then we consider the optical depth from the point 
at z{> 0) on the vertical axis {R = 0) to the upper surface 
of the disc is 
/• 00 

Ti{z) = / /9d(0)fc^exp(-/i//id)d/i = THa,oexp(-z//id), (22) 



4.4 X 10^* ergs s-^ kpc ^ ' 

Assuming the dust-to-gas ratio to be the same as in 
the Galactic disc, i.e., E{B-V) = 1.7 x lO"^^iVH(^/^0) 
(Bohlin 1978) and a Galactic extinction curve (Fitz- 
patrick 1999), we obtain the effective optical depth to 
Ha through the disc perpendicularly as 

3.3 X 1039 ergs s"! kpc'V 

Once we have r, it is straightforward to calculate the 
extinction at Ha using Equation [25] as following 



THafi = 0.5 



(28) 



At 



-2.5 log 



L 



Ha.T 



dy J ada exp(— a){exp[ — ^"''^ exp{—a/'q)y^^'^] + 
exp 



[-^exp(-a/^)(2 -//'')]}/ 



cos I 

2[l-exp(-0-$exp(-C)l, 

Similarly, replacing THa.o with rH/3,0, one obtains the ex- 
tinction at H/3. Thus Ha/H/3 ratio can be calculated, and 
then Lhq can be corrected with the extinction estimated 
from the Balmer decrement, as the extinction correction 
we make to the observed Lua- 

In order to examine if such a model can reproduce 
the observed correlations, we compute the average sur- 
face brightness Bua and color excess E{B — V) esti- 
mated by Ha/H/3 from the model based on the follow- 
ing parameters: Ehq,o, Z, cosi, rj and k. To facilitate 
comparison with observation, we obtain Bua following 
the same method used for observation (Equation [S| by 
dividing the attenuation-corrected Lh^ with the surface 
within the half-light radius, defined as the radius within 
which the Ha luminosity is half of the total Lhq. The 
E{B — V) is calculated from Ana and Anp assuming 
Fitzpatrick's extinction curve. We first fix 77 to be 1. 
Since increasing metallicity and inclination both cause 
an increase in extinction, these two factors are coupled 
in the model. Therefore, we incorporate them into one 
parameter Zi = Z/ cosi. We compute the models for 
a grid of parameters with 38.5 < log T,Ha,o < 42.0, 
Zi = 0.5, 1.5, 2.5, 3.5, 4.5, k = 0.5, 1, 2, 5.' 

The models are over-plotted on the observation sam- 
ple on the Bua versus E{B — V) plane (Figure [8]). For 
clarity, we only show the observation sample in five bins 
of Zi with a bin-width of 0.2, and on each panel the sub- 
sample is binned in Bua and shown in squares with error 
bars. We can see the observation trend on each panel, that 
E{B — V) increases as Bua increases. The model curves 
are calculated with the corresponding value of median 
Zi for each subsample. On each panel, E{B — V) gets 
larger when re increases, and each curve with constant re 
first increases and then becomes flat as Bua increases. 
At small re, i.e. the dust layer is much thinner than the 
stellar disc, E{B — V) saturates at a small value because 
the Balmer-line emission from the disc of larger height 



(29) 
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Figure 8. Average surface brightness Bjja (ergs s^^ Itpc^^, corrected for extinction) versus E{B — V), with the relative char- 
acteristic lengths of dust to Hll regions distribution in radial direction 77 = 1 (left panel) and rj = 1.2 (right panel). Observation 
sample with Z/ cosi, marked as Zi for short, in bins of 0.5 ± 0.1, 1.5 ± 0.1, 2.5 ± 0.1, 3.5 ± 0.1, and 4.5 ± 0.1 are shown. Dotted 
(blue) line represents for the relative characteristic lengths of dust to Hll regions distribution in vertical direction ft = 0.5, solid 
(red) line for k = 1, dashed (magenta) line for k = 2, and dash-dotted (green) line for k = 5. 



tends to dominate, resulting in a small E{B — V), while 
at large k more Balmer-line emission comes from the in- 
ner region close to the midplane of the disc, thus tend to 
have a large E{B — V). The saturation is due to the fact 
that the observed hght from outside weights more as the 
optical depth increases. 

Generally the k — 1 model reproduces the observed 
correlations at low Zi. However, as Zi increases, the con- 
stant K model is lower than the observation trend at the 
high-end of Bua- At medium Zi (1.5~2.5), the observa- 
tion can be better presented by the model with larger 
K (e.g. K ^ 2). This could be explained by the vertical 
expansion of dust distribution in the disc while higher 
surface brightness indicates intenser star-formation ac- 
tivity. The dust may be driven away from the midplane 
of the disc by radiation pressure or the stellar winds, 
thus the dust layer become thicker than that for low 
Bua- But E{B — V) saturates at about 0.3 mag when 
K ^ 2, which means that the effective optical depth sat- 
urates when the height of dust layer gets to twice height 
of H II regions. Another possible explanation is that the 
Kennicutt-Schimidt law we adopt in the model is an av- 
eraged star formation law that linked between the sur- 
face density of gas and star formation as the form of 
SsFR oc E^. Leroy et al. (2008) suggested that the star 
formation law may vary as a function of local conditions 
at different galactocentric radius, including metallicity 
and the dominated gas content (Hi or H2). The SF law 
has smaller A'^ value at higher Esfr in the inner region 



of a galaxy (e.g. Bigiel et al. 2008, their Figure 11), thus 
large dust reddening in the inner region dominates the 
observed E(B — V\ As a result, the dust reddening is 
higher than that the present model predicts. 

In the model above, we assume 77 = 1, i.e. the dust 
and Hll region has the same characteristic scale in the 
radial distribution, but the gas/dust may distribute at 
larger radii than the star-forming regions (e.g. Bigiel et 
al. 2008). Thus we also consider 77 = 1.2, and find that 
the observation with large Zi can be well explained by 
model with k = 1 at low Bhq and k ^ 2 at high Bna- 
It is possible because in a spiral galaxy centre, where the 
metallicity is higher than outside, we are observing the 
star-formation dominated centre, and the gas (or dust) 
distribution is more extended. 

In this simple toy model, the optical depth rua is as- 
sumed to be proportional to the metallicity Z. It comes 
from the assumption that pd oc pg x Z. However, there 
is evidence that this assumption may be not the case in 
practice. Boissier et al. (2004) studied six nearby late- 
type galaxies in FIR and UV images, and find that the 
dust-to-gas ratio is proportional to (^/Z©)"'**, which was 
fiatter than the linear correlation we adopted above. This 
suggests that the dust-to-gas ratio may not evolve lin- 
early with Z. In addition, we assume in the model that 
the disc height is very thin, i.e. the disc height is quite 
small compared to the disc radius, thus the radial gra- 
dient in gas density can be neglected when calculating 
THct- Under that assumption, the effect of inclination on 
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Figure 9. Same figure as Figure [S] except tiiat the model is 
calculated with Equation 1281 and 1291 but adjusting Z/Zq to 
{Z/Zq)°-'', and cosi to (cosi)"''^. 



optical depth can be expressed as (cosi)~^. However, in 
practice, the disc may not have a uniform height, for in- 
stance, it may be thicker in the centre and thinner at 
large radii, or the presence of bulge in the disc centre. 
Also, the distribution of dust and star-formation regions 
may be not smooth, for example when the spiral arms 
are present in the disc galaxy. In those cases, the incli- 
nation factor might be flatter than (cosi)"'^. Accounting 
for these realistic effects, we adjust Z/Zq to (Z/Zq)" in 
Equation l28l and (cosi)~^ to (cosi)~" in Eauationl25Jand 
1291 then test several values for the power-law indej{j. It 
turns out that with a = 0.7 the simple model reproduces 
the observation trend well (Figure [9]). When Zi is low, 
the model of constant k (0.5 or 1) reproduces the obser- 
vation trend well. While Zi is higher, the observation can 
be reproduced with the model of k ^ 2. 

To summarize, the simple toy models reproduce the 
observed trends between dust reddening and Ha surface 
brightness, that reddening increases with Sua- The com- 
parisons between observation and the toy models indicate 
that, the relative distribution of Hll regions to gas/dust 
both in vertical and radial direction of the galaxy disc, 
may vary in different environment of star-formation ac- 
tivity levels or metallicities. The dust might be dis- 
tributed further beyond the Hll regions when the star 
formation activity is intenser, or the dust is distributed 
in a thicker disc (k ^ 1) than the Hll regions when the 



^ Because the metallicity and inclination are coupled in this 
simply toy model, the power-law index of cos i is adopted to 
be the same as that of Z, so as to be incorporated into Zi. 



metallicity gets higher. In addition, our results imply that 
dust-to-gas density ratio may be not proportional to the 
metallicity Z, but to a power-law of Z, Z" with a --^ 0.7. 



6 CONCLUSIONS AND IMPLICATIONS 

We present the empirical formulae of intrinsic reddening 
as a function of intrinsic Ho? luminosity or surface bright- 
ness, metallicity and disc inclination fEguation \7\ 1 1 II and 
Table [2] and |3]) for a large sample of '~-^22 000 well-defined 
star-forming disc galaxies selected from the SDSS. With 
the empirical formulae, the reddening parameterized by 
E{B—V) could be predicted within la uncertainty of 0.07 
mag. The AE{B — V), defined as the observed reddening 
estimated from the Balmer decrement minus predicted 
value, does not correlate with the three parameters used 
in the formulae. We also find that AE{B — V) is inde- 
pendent of the stellar mass, 4000 A break strength or 
electron density within the galaxy. 

The observed trend between the dust reddening and 
Ha surface brightness could be reproduced by a plane- 
parallel slab toy model, in which the dust is scaled with 
gas, and all the distributions of the dust, gas and star- 
formation regions are smooth and follow exponential laws 
in the disc. We find from the comparisons of the models 
to the observation, the relative vertical scale of dust dis- 
tribution to Hll regions distribution may vary with Ha 
surface brightness or metallicity. The higher intensity of 
star formation could drive the dust to distribute in a 
thicker layer, or the dust disc is thicker than the disc of 
the H II regions when the metallicity is higher. However, 
the real galaxies may be not fully described by this sim- 
ple toy model. For example, when the spiral arms are 
present, or the SF regions are clumpy, the assumptions 
of smooth distributions will be improper. In addition, 
our result implies a different scaling law of dust-to-gas 
ratio as a function of metallicity Z from the linearly re- 
lation (Equation [TSj. Our empirical relations, which sug- 
gest dust reddening is partially dependent on metallic- 
ity (cx ^0-6~o.7^^ observational constraints on 
the output of radiative-transfer model calculations, and 
hence provides constraints on the assumed dust-to-gas 
ratio in the model. 

Since metallicity-dependence has been introduced, 
we believe that our empirical formulae should be appli- 
cable to the early evolutionary stage of a galaxy, when 
the gas metallicity might be low. In the starburst galaxies 
at high redshift, dust surrounding star-forming region has 
been proven to be prone to sub-milimeter and far-infrared 
observations. But most of those galaxies are relatively 
dim in optical, implying severe extinction and reddening 
to these galaxies. The observed galaxies at high redshifts 
are affected by reddening, in the way that the galaxies 
with heavy attenuation and reddening may escape from 
the detection. The importance of the reddening correc- 
tion depends on the wavelength band that is used. For 
example, UV continuum luminosity is usually used to in- 
fer SFR in high-redshift galaxies. It is very sensitive to the 
attenuation/reddening. Calzetti et al. (2000) find that the 
extinction to emission lines and to stellar continuum in 
optical band are different, and the continuum reddening 
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is only about 40 percent of that for emission lines. How- 
ever, the extinction to UV continuum is likely similar to 
the H II region because UV continuum is also emitted by 
young hot stars, which are likely surrounded by Hll re- 
gions. As a result, our results may be used to estimate 
the reddening for these high-redshift galaxies with cau- 
tion, that the metallicity should be well-calibrated and be 
within the range covered by our sample. Following Equa- 
tion [7] and the discussions in Section |4l SFR in luminous 
star formation or high metallicity galaxies might have 
been under-estimated (e.g. Heckman et al. 1998; Hop- 
kins et al. 2001; Panuzzo et al. 2007), or those galaxies 
are even completely lost in optical surveys. 

The quantitative relationship of reddening with the 
metallicity, luminosity or surface brightness, and the in- 
clination for disc galaxies can also be incorporated into 
the current semi-analytic models of galaxy formation and 
evolution (e.g. Croton et al. 2006; De Lucia et al. 2004; 
Kang et al. 2005). In those models, a simple chemical 
enrichment scheme has been included, and the metal 
abundance has been actually predicted. With the known 
metallicity and the star formation rate (intrinsic lumi- 
nosity) provided by the models, the reddening can be 
estimated accurately with our formulae, and then used 
to shape the spectral energy distribution and predict the 
emerging luminosity more realistically. 
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